function C = C_plane_stress(E,nu)
C = E/(1-nu*nu)*...
    [   1 nu 0;
        nu 1 0;
        0 0 (1-nu)/2    ];
end
